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We give an analytic demonstration that the 3 + 1 dimensional large SU{N) pure 
Yang-Mills theory, compactified on a small so that the coupling constant at the com- 
pactification scale is very small, has a first order deconfinement transition as a function of 
temperature. We do this by explicitly computing the relevant terms in the canonical par- 
tition function up to 3-loop order; this is necessary because the leading (1-loop) result for 
the phase transition is precisely on the borderline between a first order and a second order 
transition. Since numerical work strongly suggests that the infinite volume large A^ theory 
also has a first order deconfinement transition, we conjecture that the phase structure is 
independent of the size of the S^. To deal with divergences in our calculations, we are led 
to introduce a novel method of regularization useful for nonabelian gauge theory on S^. 
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1. Introduction 

It is widely believed that 3+1-dimensional SU{N) Yang Mills theory on confines at 
low temperatures, but is deconfined at high temperatures. Compelling numerical evidence 
indicates that in the absence of quarks, when all fields are in the adjoint representation, 
there is a sharp phase transition separating the confined and the deconfined phases, which 
occurs at a temperature T ~ Aqcd- Since the Yang Mills theory is strongly coupled at the 
transition temperature, the deconfinement phase transition is rather poorly understood. In 
particular, using the currently available analytic techniques it is not possible to determine 
even the order of the transition for 7^ 3; lattice simulations suggest that the transition 
is of second order for N = 2 and of first order for > 3 (see [^|2| for the latest results for 
N > 3). 

The intractability of the thermal behaviour of Yang Mills theory on IR"^ is related to 
the absence of a dimensionless coupling constant. It is thus interesting to note that Yang 
Mills theory compactified on an of radius R has an effective dimensionless coupling 
constant, RAqcd- Indeed, when RKqcd ^ 1, the Yang Mills coupling constant is weak 
even at the lowest energy scale in the theory, i? ~ As a consequence, at small values 
of AqcdRj the thermal behaviour of this system is completely tractableB. Unfortunately, 
the most interesting feature of infinite volume thermodynamics - the sharp deconfinement 
phase transition - is smoothed out into crossover behaviour at any finite R, assuming that 
is also finite. 

However, in the 'thermodynamic' ^ 00 't Hooft limit (with fixed A = 9ym^ 
0), this deconfinement phase transition remains sharp even at finite R. In this limit it is 
possible to study the dynamics of the deconfinement transition as a function of the effective 
coupling constant RKqcd- When RAqcd ^ 1, this system approaches the theory on IR"^. 
On the other hand, in the opposite limit RAqcd the theory is weakly coupled, and 
may be solved exactly quite remarkably it turns out that this 'free' gauge theory 

undergoes a confinement- deconfinement phase transition at a temperature of order 1 / 

^ More precisely, it is tractable for temperatures smaller than an upper bound that scales to 
infinity in the limit AqcdR 0, see section 5. This would not be true if the gauge field had zero 
modes on the compact space, which is why we chose a sphere rather than, say, a torus. 

^ As discussed in the A = theory must still obey a Gauss Law constraint which requires 
physical states to be gauge invariant. This constraint leads to nontrivial thermodynamics even at 
zero coupling. 
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At strictly zero coupling, the transition is first order, but lies precisely at the border 
between first and second order behaviour, as reviewed below. Consequently, to understand 
the nature of the transition at weak nonzero coupling, the leading effects of interaction 
terms must be taken into account via a perturbative calculation. This calculation is the 
goal of the present paper. 

Before describing the calculation and our result, we recall the essential details of the 
story in the A = limit. It was demonstrated in [^,|^ that in the limit RAqcd ^ in 
which the theory becomes free, the thermal partition function of Yang Mills theory on a 
3-sphere of radius R reduces (up to an overall constant) to an integral over a single unitary 
SU{N) matriJ 

Z{P)= j[dU]eM-SeffiU)], (1.1) 



where 



Sef^iU) = - f; '^^'"J'^''^ TV(^-)IV(^-); zvix) = (1.2) 



n=l 



and (3 = 1/T. The matrix U is the holonomy of the gauge field around the thermal circle 
(more precisely, U = e*^" with a the zero mode of Aq on S'^ x [^). All other modes of 
the theory are massive, and the effective action is obtained by integrating them out. 

As usual, in the large limit, it is convenient to replace the integral over the unitary 
matrix U by an integral over the eigenvalue distribution p{9) = J2i ^{^ — ^i) (where e*^', 
for i = 1, ■ ■ ■ , A, are the eigenvalues of U). Let Un denote the n^^ Fourier mode of the 
eigenvalue distribution, Un = J e^'^^ p{9)d9 = Tr(t/")/A. In the large A limit, ( |1.1|) may 
be rewritten as IHl 



Z{l3) = / dundun exp 



(l-;zy(e-"^/^)) 



Ur. 



1 

n=l 



:i.3) 



with a complicated integration boundary for the u^s coming from the non-negativity of 
p{0). Equation ( pT3|) may be evaluated in the saddle point approximation in the large A 
limit. This system has one obvious saddle point located at w„ = (for all n > 1). This 
saddle point is stable and dominates ( |1.3| ) whenever (1 — zv{e~'^^^^)) is positive for all n, 
which is the case for T <T^= ^ a7M326_ 



We will generally ignore the distinction between SU (N) and U (N) groups in this paper; the 
only difference between their partition functions is an overall coupling-independent factor coming 
from the free U{1) photons. 
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At T = Tc, zv{e~^^^) = 1 and the potential for ui in (L3) is exactly flat (all other 
Un remain massive). For T > Tc, the variable ui becomes tachyonic about the saddle 
point described above. At these temperatures, the system is dominated by a new saddle 
point, one in which ui has an expectation value of order unity (see 0). Thus, free Yang 
Mills theory on undergoes a phase transition at T = Tc; the order parameter for this 
transition is the expectation value of the Polyakov loop iti (more precisely, since ui has an 
arbitrary phase, the order parameter is actually [H). Since the saddle point changes 
discontinuously at T = T^, this phase transition is of first order. However, this phase 
transition is extremely finely tunedS, in a sense we will now explain. 

As described above, the exact Wilsonian effective action for the order parameter ui, 
expanded about the low temperature saddle point in the limit of zero gauge coupling, 
is Seff = N"^ (l — -2\/(e~^/^)) This effective action is corrected at nonzero gauge 

coupling; it was argued in [|| that at nonzero coupling it takes the formi 



S, 



eff 



Ar2 



(l-zv.(e-^/«)) +0(A)] |tii|2 + A26(/3)|tii|4 + 0(A4), (1.4) 



where A = Qym^ Hooft coupling evaluated at the energy scale l/R, and b{(3) is 

a perturbatively computable function of the temperature. 

The action ( [1.4| ) describes a system that undergoes a phase transition at T = Tc+0{X). 
The nature of this phase transition depends crucially on the sign of 6 = b{Pc), where 
Pc = ^/Tc- If h is negative, the phase transition is of first order, as in the free theory. 
However, unlike in the free theory, this transition occurs at a lower temperature than the 
temperature at which ui in ( |1.4| ) becomes massless (the latter temperature was identified 
in with the Hagedorn temperature of the large Yang- Mills theory). 

On the other hand, if h is positive, the phase transition continues to occur precisely 
at the temperature at which ui becomes massless; however, it is now of second order and 
is followed, at a slightly higher temperature, by another phase transition of third order, 
similar to that of (see section 6 of for more details). 

The leading perturbative contribution to the value of h at the phase transition tem- 
perature is determined by a set of 2-loop and 3-loop vacuum diagrams in the Yang Mills 



^ As one indication, the order parameter ui is massless at the phase transition point, a feature 
usually associated with second order transitions. 

^ To obtain the explicit effective action, we first integrate out all massive modes on the sphere, 
and subsequently all of the other modes for n > 1. 
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theory on a sphere 0. In this paper we evaluate b by computing the relevant Feynman 
diagrams. 

Our main result is that b{Pc) — —5.7 x 10~^. Note, in particular, that b is negative; 
consequently, the deconfinement transition for large SU{N) Yang-Mills theory on a 3- 
sphere with small but nonzero RAqcd is of first order. As noted above, lattice simulations 
suggest that the large deconfinement phase transition is also of first order in the opposite 
infinite- volume limit RAqqd oo. Thus, it is tempting to conjecture that the phase 
diagram for the large SU{N) Yang Mills theory on takes the form shown in figure 1. 



Hagedorn Temperature 



TR 




CONFINED PHASE 



RAoc 



Figure 1: The simplest possible phase diagram for large Yang Mills 
theory on S^. 



The solid line in figure 1 denotes the phase transition. The dotted line describes 
the boundary of stability (the locus at which the coefficient of is zero) of the low 
temperature phase, which would be interpreted by a low temperature observer as a limiting 
or Hagedorn temperature. 

Of course, the conjectured phase diagram in figure 1 merely represents the simplest 
phase diagram consistent with our knowledge of the behaviour of all order parameters at 
weak and strong coupling. It is possible that the true phase diagram is more complicated; 
for instance, the confined phase at small RAqcd could be separated from the confined 
phase at large RAqcd by a phase transition. 

This paper is organized as follows. In section 2, we set up our calculation, and describe, 
in general terms, the procedure we employ to compute b. We then proceed in section 3 to 
enumerate and evaluate all diagrams that contribute to b. 
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The lengthy calculation we present is straightforward in principle but involves some 
subtleties. First, it is necessary to deal with UV divergences which show up in the perturba- 
tive computation. Due to technical difficulties associated with implementing dimensional 
regularization on S^, we have found it convenient to regularize the Feynman diagrams in 
a rather unusual fashion: we use a (non-gauge-invariant) momentum cut-off, and simul- 
taneously add in a set of compensating non-gauge-invariant counterterms, such that the 
full theory is gauge-invariant when the cutoff is removed. This regularization procedure is 
described briefly in section 2.4, with a more detailed discussion, including checks, exam- 
ples, and the explicit evaluation of the counterterms needed for our calculation presented 
in section 4 and in appendix A. 

In section 5, we discuss potential problems related to the infra-red divergences asso- 
ciated with finite-temperature field theory. Naively there should be no trouble at finite 
volume, where all modes are massive. However, it turns out that there is still a breakdown 
of perturbation theory at sufficiently high temperatures, when the dynamically generated 
mass scales exceed the Kaluza-Klein scale 1/R. Fortunately, as we argue in section 5, these 
effects are not important at the transition temperature and do not affect our determination 
of b. 

In section 6, we present some conclusions and discussion. 

Finally, in order to compute sums over spherical harmonics, which replace the loop 
integrals of the flat space theory, we needed to derive various spherical harmonic identities. 
These identities, together with some basic properties of the spherical harmonics for are 
collected in appendix B. 

2. The Setup for the Perturbative Calculation 

As we have described, the order of the deconflnement phase transition for pure Yang- 
Mills theory at small volume is determined by the sign of the quartic coefficient b in the 
effective action ( p..4| ) at the A = transition temperature. In this section, we set up the 
calculation of this coefficient, which may be determined from leading order perturbative 
corrections to the matrix model action ( |1.2|) for the SU{N) pure Yang-Mills theory on 
S^. For simplicity, we will generally take the radius of the to be one; it can always 
be reinstated by dimensional analysis. The actual computation is presented in the next 
section. 



5 



2.1. Basic objective 

The basic set-up for the computation was described in section 4 of 0. We consider 
pure SU{N) Yang-Mills theory on S^, at finite temperature. The thermal partition func- 
tion may be computed by evaluating the Euclidean path integral with Euclidean time 
compactified on a circle of radius (3 = 1/T. The Euclidean action is given as usual by 
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/:=-y dt J d''xtr{F^,F^"'). (2.1) 

For calculations on S^, it is convenient to work in the gauge 

d^A' = 0, (2.2) 

where i = 1, 2, 3 runs over the sphere coordinates, and di are (space-time) covariant deriva- 
tives. Equation (|2.2| ) fixes the gauge only partially; it leaves spatially independent (but 
time dependent) gauge transformations unfixed. We fix this residual gauge invariance by 
setting the constant mode of Aq to be constant in time, 



dt Ao = 0. (2.3) 

It will be convenient to give this (time independent) zero mode a name; we define 

a = ^ / Ao, (2.4) 

where ^3 is the volume of the 3-sphere. a will play a special role in what follows, because 
it turns out that it is the only zero mode (mode whose action vanishes at quadratic order) 
in the decomposition of Yang Mills theory into Kaluza-Klein modes on x S^. 

As a is a zero mode, it cannot be integrated out in perturbation theory (roughly 
speaking, a fiuctuations are always strongly coupled in the bare action). In order to 
perturbatively evaluate the free energy we will therefore adopt a two step procedure. In 
the first step we integrate out all nonzero modes and generate an effective action for a. As 
described in , this action will be non-trivial even at zero coupling, and it is corrected in 
perturbation theory in A. In the second step, we analyze the remaining integral over a. 

On general grounds described in 0, the finite temperature effective action for a can 
be written completely in terms of the unitary matrix U = e*^" in the form : 

Seff=Yl Cm,-m(x)tr([/-)tr(t/-"^) + A/3 Cm,n,-n.-nixMU"')triU^)triU-^-^)/N 

m m,n 

+ A2/3 C^,n,p,-m-n-p(x)tr(Otr(t/")tr(t/^)tr(t/-"^-"-P)/iV2 + ---, 

m,n,p 

(2.5) 
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where x = e~^. Here, the coefficients obey appropriate constraints such that the action is 
real, and in general are corrected at higher orders in A and As described in [Q, the 

free energy F of the Yang Mills theory is then given by the matrix integral 

e-f^^ = [ [dU]e-^'^^f^^\ (2.6) 



where the Fadeev-Popov determinant corresponding to the gauge-fixing (^]^) transforms 
the integral over a into an integral over the gauge group with Haar measure [dU]. 

At large N, the unitary integral may be evaluated using saddle point techniques. 
Defining Un = tr(t/")/A^, we write the effective action in the form 

[dw,]e-^'^^//(^"). (2.7) 

where S'^j-j- includes contributions both from Sef / and from the Vandermonde determinant 
obtained in changing to the variables Un- The order A^^ contribution to the free energy 
is then given by the minimum value of S'^j-jr{un), and the deconfinement phase transition 
occurs where this minimum is no longer at \un\ = 0. 

As described in detail in section 6 of 0, in order to compute the order of this phase 
transition, we have to look at S'^^j{un) near the phase transition point Xc = 2 — where 
(as shown in and reviewed in the introduction) the mass term of ui changes sign, and 
compute the leading corrections to the potential for ui. The relevant terms in the action, 
to leading order in A and in x — Xc, take the form 

S'es{Un) =/il(Xc - x)|wi|^ + ^12\U2\^ + • • • 

+ A/3[Ci,i,_2(^x?^2 + u\u2) + ■■■] (2.8) 
+ A2/?[Ci,i,_i,_i|^ii|4 + ■■■]. 



At large A, the effective action ( |1.4| ) for ui is obtained from this by classically minimizing 
over Un for fixed ui. In particular, the variable U2 may now be integrated out in (|2.8| ) by 
setting it to its classical value 

„2 ^ -aMIzIz^^s ^ 0(A2). (2.9) 
This yields the following effective action for ui : 

S',fj{ui) = ^lixc - x)|wi|2 + 6AVi|^ + ■ ■ ■ , (2.10) 



where ^ 2 

b = /3cCi,i,_i,_i - ^2£hh^. (2.11) 

Our goal will be to compute the coefficients appearing here at the leading order in A, and 
thus determine the sign of b, and the order of the transition, at weak coupling. Since an 
71-loop diagram has at most n + 1 index lines, a term in Seff with m traces gets its lowest 
order contributions at m — 1 loops. Thus, the term fX2 arises at one-loop order, Ci,i,_2 
requires a two-loop computation, while Ci^i^_i^_i requires a three- loop calculation. 

2.2. Gauge-fixed action 

We now set up the perturbative computation that will determine Sef / (U) , the effective 
action for U (which we treat as a background field). The Fadeev- Popov determinant 
corresponding to the gauge fixing condition (|2.2|) is 



detd,D' = J VcVce-''^^^^^'^\ (2.12) 

where denotes a gauge covariant derivative 

DiC = di- igYM[Ai, c], (2.13) 

and c and c are complex ghosts in the adjoint representation of the gauge group. The 
quadratic terms in the gauge-fixed Yang-Mills action ( |2.1| ) take the form 



- j d^xii (y^MDl + d'')A' + ^Aoa^Ao -f ca^c^ , (2.14) 

where is the Laplacian on the sphere and 

DrX = doX -i[a,X]. (2.15) 
The interaction terms in ( p.l| ) are given by 

j d^xti{igYMDrA'[Ai, Aq] - igYM[A\ A^]diAo - igYMdiAj[A\ A^] + 



- ^[Ao.A,][A^,A^]-igYMd.c[A.c]). 



2.16) 



2.3. The spherical harmonic expansion 

On S^, integrals over spatial momenta are replaced by sums over the quantum numbers 
of SO [4) angular momenta. It will thus be useful to expand the fields explicitly in terms 
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of an orthonormal basis of functions on which are angular momentum eigenstates, 
and write the action exphcitly in terms of a standard set of integrals over these spherical 
harmonic functions. We denote the scalar and vector spherical harmonics on by S"{Q) 
and Vf{fl), where a = {ja, nia, m'^) and (3 = [jp, m^, m'^, ep) are the angular momentum 
(and parity for the vector field) quantum numbers for the various modes. The properties 
of these functions are reviewed in appendix B. 

We expand the modes of the fields in terms of these spherical harmonics as follows : 

a 

A,{t,e) = Y^AP{t)vt{e); (2.17) 
c{t,e) = Y,c''{t)s%e). 

a 

Note that general vector functions also include modes proportional to VS"", but these are 
eliminated by our gauge choice. Below, it will be useful to denote the complex conjugates 
of and by S"^ and V[ . 

In terms of these spherical harmonics, we define 

jj^Pi ^ /" . yPs^, (2.18) 
is3 



I S3 

where explicit expressions for C, D and E may be found in ^ and are collected in appendix 
B. Note that C is antisymmetric in a and 7, D is symmetric in a and (3, and E is totally 
antisymmetric. 

Using the spherical harmonic expansions, we may now write the action for gauge- fixed 
pure Yang-Mills theory on explicitly in terms of modes. The quadratic action becomes 

S2 = j dttr(^A^{-Dl + (i„ + lf)A" + ^a^'j^U^ + 2)a" + c^j^{j^ + 2)c"). (2.19) 

In addition, we have cubic interactions 



S3 = gvM j c^nr(zc"[A'^, c^]C"^^ + 2ia"A'^a^C"^^ 



(2.20) 

D^A^ja^D"^'^ + iA'^A^A^e^ii^ -f 1)£;"^^), 
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and quartic interactions 

J ^ V J\[Jx + ^) J (2.21) 

The propagators of the various fields follow from ( |2.19| ) and are given by 

(C(i)cf,(t')) = ■ ./^./ "^^(^ - t')SadSct, (2.22) 

= ■ ^/^o/ "^^(^ - t')^ad5ci„ (2.23) 

Ja\Ja "r 

(^a.(t)^f.(t')) = 5"^A,t'^'(t - t')- (2.24) 
Here, A is defined to be a periodic function of time satisfying 

{-Dl + {j + lf)/\,{t) = 5{t) (2.25) 

where we have suppressed matrix indices. For < t < /9, the explicit solution is given by 



with the value for other values of t defined by the periodicity. Here, a is shorthand for 
a (g> 1 — 1 ® «, and a term a"^ (g) a'^ in the expansion of A should be understood to carry 
indices (a")"^(a"")^^ in {^^. 



For our calculations, the following correlators are also useful 

{DrA:,{t)A^^,{t')) = -{A^.mX.it')) = 5-~^D^/\ff\t - t% (2.27) 
{D^A':,it)D^A^^^{t')) = r^5(t - t')6adScb - r^(j« + lfAff\t - t'). (2.28) 

2.4- Regularization and counterterms 

As usual with four dimensional gauge theories, certain perturbative calculations lead 
naively to ultraviolet divergences. In our calculation of the coefficient b in ( |1.4D , we will 
find that all divergences cancel, but only after summing a collection of logarithmically 
divergent diagrams. It is therefore necessary to introduce a regularization scheme, and in 
order to obtain the correct finite result, this must respect gauge invariance. 
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In principle, there is no obstacle in taking the usual approach of applying dimensional 
regularization. In our case, this amounts to considering gauge theory on 5""^ x S"^ x IR'^~^.il 
While it is straightforward to write down the appropriate expressions for Feynman dia- 
grams in this dimensionally regulated theory, we have found it difficult to evaluate them 
because of the combination of momentum sums and integrals that appear. Thus, we have 
found it helpful to apply a more unusual approach, which nevertheless gives precisely the 
answers that would have been obtained through dimensional regularization. 

In practice, we apply a sharp momentum cutoff to the total angular momentum quan- 
tum number for modes on the sphere. This does not respect gauge invariance. However, 
as explained originally by 't Hooft Q in his series of classic papers on the renormaliz- 
ability of Yang Mills theory, a non gauge invariant regularization yields gauge invariant 
results when employed with a bare Yang Mills action that includes an appropriate set of 
non-gauge-invariant counterterms. As we describe in section 4, these counterterms may 
be determined by demanding that simple Green's functions evaluated using the cutoff and 
counterterms agree with the same Green's functions evaluated using dimensional regular- 
ization. Apart from curvature-dependent counterterms (which do not contribute to our 
calculation) , this comparison may be carried out in flat space, since all counterterms must 
be local. 

In section 4, we present the calculation to determine the precise coefficients for all 
counterterms necessary in our calculation, together with a more complete discussion of the 
regularization scheme and a variety of consistency checks. In the end, we should expect 
non-zero coefficients for all counterterms with dimension less than or equal to four which 
respect SO{3) invariance. 



3. The Perturbative Computation 

In this section, we proceed to calculate the coefficients in ( |2.8[ ) necessary to determine 
the order of the deconfinement transition at weak coupling. We have 

J (3.1) 

q1 loop „ „ 

^ As we describe in section 4, in order that all momentum sums/integrals are rendered finite 
when employing dimensional regularization in Coulomb gauge, it is necessary to analytically 
continue both the number of dimensions which participate in the Coulomb gauge condition, and 
the number of dimensions which do not, a procedure referred to as split dimensional regularization. 
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where the expectation value in the last line is evaluated in the free theory with propagators 
given in section 2. The required leading order contributions to /i2, Ci,i,_2, and Ci,i,_i,_i 
appear at one, two, and three loops respectively. 



3.1. Simplifying the action by integrating out Aq and c 

Since the action is quadratic in a and c, these may be integrated out explicitly to yield 
additional interaction vertices for the A's. The first contribution arises from loops of a or 
c. For the calculation up to three loops that we are interested in here, the relevant vertices 
in the resulting effective action (combining a and c loops) are a quadratic vertex 

J/3U/3 + 2) 



a cubic vertex 



f^a-yiP T-)-f2Sf3 T-j-ysSa 

As = -2igl^N5{Q) . M A^-A^-A^% (3.3) 

Ja(ja + 2)j^(j^ + 2) 

and a quartic vertex 

((J(^liP(Jl3l2P ]J13>^P 2)74^(5 ^ 2^72A/3 jj-y^pP ]Jl4pa 

^Uj» + 2)jp{jp + 2)jp{jp + 2) + 2 j^(j^ + 2)jf3{jp + 2) 

X (^tr{A^'A^^A^^A'^^) + j^tr{A'^' A^^)triA'^^ A'^^) 

+ -^tr(A^M^3)tr(A^M^4) + -^tr(A^M^4)tr(A^2A^3) 

(3.4) 

We have included here only terms that can contribute to planar diagrams. Note that all 
of these are proportional to a divergent factor 5(0), and so any diagrams containing these 
vertices must eventually cancel.0 Note also that since these vertices arise from loops, they 
have additional factors of gym-^ compared to ( p.20| ), (|2.21| ). Therefore, inserting any one 
of these vertices into a diagram counts as an additional loop. 

In addition to these, we have vertices arising from open strings of a's containing two 
vertices linear in a and some number of vertices quadratic in a. These start at quartic 



^ The divergence associated with 5(0) terms may be regulated using a momentum cutoff in the 
direction. We have checked that in the properly regulated theory, the naive cancellations of 
5(0) terms described below persist without introducing any new finite contributions. 
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order, and for our three loop calculation, we will need the quartic, quintic, and sextic 
vertices. These are 

SyM ( _D"i/^l'^(7*7iT(^T72A£^a2/32A JJotiPia jj-y-^Xa Jj'y2>^f jja2f32T 

" ^ \ jA(iA + 2)j,(i, + 2)jV(jV + 2) ^ iaO. + 2)jVOV + 2) 
tr([[A"S D^A^i], I^^A^^]^ ^72])_ 

Using these effective vertices, it is straightforward to check that all divergent contributions 
proportional to 5(0) cancel. Any diagram with a vertex A2 will cancel the 5(0) part of 
a similar diagram with A2 replaced with a B4 having its two covariant derivative legs 
contracted. Similarly, any diagram with a vertex A3 will cancel the 5(0) part of a similar 
diagram with As replaced with a S5 having its two covariant derivative legs contracted. 
Finally, any diagram containing A4 will cancel against a combination of two diagrams: one 
with A4 replaced by Bq with its two covariant derivative legs contracted, and one with A4 
replaced by two B4S with the four covariant derivative legs contracted into a loop. 

Thus, all 5(0) terms coming from the An vertices cancel out, and it is easy to check 
that there are no additional 5(0) terms coming from diagrams with B^ vertices apart from 
those needed to cancel the A^ vertex diagrams. As a result, we may proceed with the 
calculation by keeping only the transverse photons A", and evaluating all diagrams built 
from the original cubic and quartic vertices in the transverse photons plus the additional 
vertices S4, S5, and i?6, ignoring any terms proportional to 5(0). 

The diagrams contributing to the free energy at one, two and three loop orders, after 
having integrated out Aq and c, are shown in figure 2. The S-type vertices are denoted 
by circles. 

3.2. 1-Loop 

The one-loop computation of the path integral was described in § . The result (writing 
only the leading terms in the large limit) is 

Seff{U)= — - — tr{U )tr{U ), (3.6) 

n=l 
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One loop: ^^^^^^^^^ 
la 



Two loops 




Three loops 





Figure 2: The diagrams contributing to the free energy up to 3-loop order. 
In this figure we present a particular planar form for each diagram, but in 
some cases the same diagram may also be drawn in the plane in different ways. 
There is also an additional counter-term diagram that will be discussed later. 

where zv{x) is the single-particle partition function for a free vector field on S^, given by 

6x^ - 2x^ 



zv{x) 



(3.7) 



(1-X)3 

Changing variables to Un = tr {U'^)/N and including the additional Vandermonde deter- 
minant from the measure, we have 



^^//K) = ^ + E ^(1 - .v^M)KP 



n=l 



and thus 



Ii2{xc) = -(1 - zvixD) ^ 0.481125 . 



(3.8) 



(3.9) 
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3.3. Two Loops 

In this section we compute the coefficient Ci,i,_2 in ( |2.(j| ) by evaluating the three 
two-loop diagrams of figure 2. Note that since the propagators depend only on the relative 
time between vertices, it is always possible to change variables so that the integrand is 
independent of one of the time variables. The integral over this variable then gives an 
overall factor of 13 so it is convenient to define F{U) = Sef f{U) / p. In terms of the 
propagators and spherical harmonic integrals defined above, we find that the three two- 
loop contributions to F are! (with summation over the spherical harmonic indices a, P 
and 7 implied) 



'2a 



i7(i7 + 2) 



F2c — 9ym ■ /■ 



(3.10) 

where we have defined 

^ E»P^[e^[j^ + 1) + ep{jp + 1) + e^(i^ + 1)) . (3.11) 

In the expressions above, each of the propagators contributes factors of a to two of the 
three index loops, which we label by a, 6, and c. The notation a^b indicates that for the 
tensor products {a®l — l®a) appearing in the propagator, the first and second elements 
of the tensor product appear in the traces associated with index loops a and h respectively. 

The expressions ( |3.10|) involve sums over products of spherical harmonic integrals. In 
all cases here and below, the sums over quantum numbers m, m' and e may be carried out 
explicitly using the formulae in appendix B. To express the results, it is useful, following 
, to define some functions which appear in integrals of products of spherical harmonics 



R2{x,y,z) 



{-ly' r(x + 1)(^ + l)(a' - x){a' - y){a' - z){a' + 1)^ 



TT 



+ 



(3.12) 



Recall that in diagram 2c, we ignore the part proportional to 5(0). 



15 



{{e^ix + 1) + e,{z + l) + y + 2)(e^(x + l)+€,{z + l)+ y) 
{e^{x + 1) + e^{z + 1) - y){e^{x + 1) + e^{z + l)-y-2))^ , 
(-1)"'+' . , , f{a' + l){a' -x){a' -y){a' - z)Y 



1 

2 



^e^eyc^y J ^ & Kx y zjy 4(x + 1) (y + 1) + 1) 

• ((e^(x + 1) + ey{y + 1) + e^iz + 1) + 2)(e^(x + 1) + ey{y + 1) + e^(z + 1) - 2)) 
R^e^^ye, {x, y, z) =RAe^eye, (x, y, z){e^{x + 1) + + 1) + + 1) ) , 

(3.13) 

where the right-hand sides of the equations are defined to be non-zero only if the triangle 
inequality \x — z\ <y< x+z holds, and if a = {x+y+z)/2 (in R^) and a' = {x+y+z+l)/2 

(in R2 and R4) are integers. We also define R^-^. = R3++ = R3--, R3- = -Ran = R3 1-, 

i?4-(- = -R4-1 — \ — |- and R4— = R4-\ — | — 

With these definitions, we find after performing the sums over m, m' and e, that 



F2a = ^|^Ja(ja + 2)j^(i^ + 2) A, JO, tt^fe) A,^ (0, , 
1 



1 

2 d2 / ■ ■ „• ^ I / „• I „• „• I 1 ^ 2 d2 



F2b = -9YMi-^ija+jp+jl + 3) ^4+(ia, J/3, J7) + ija+j/3-j-, + 1) ^4- O'a, i/3, J't)) 



Aj^ (t, ttafe) (t, ttfoc) (t, ac 



F2c = -^^^{Rl+{jc.,jjJ(3)+Rl-{jc.,jj,jf3)) 



, CK&c) Aj^ (0, Ctab))- 

(3.14) 

These expressions are all to be summed over the j's, with the sums unconstrained in F2a, 
and constrained in F2b and F2c by the rules given above. 

As described above, in order to analyze the phase transition we need to compute the 
specific term in F2a + F2b + -p2c of the form 

Ci,i,_2(7^^,(tr(t/)tr([/)tr((t/t)2)+tr(t/2)tr([/t)tr([/t)), (3.15) 

and to determine the coefficient Ci,i,_2 at the deconfinement temperature of the free 
Yang-Mills theory, Xc = 2 — v^. For each diagram, we therefore expand the product of 
propagators in powers of U, and sum the coefficients of all terms of the form (|3.15|) . 
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For F2a, all sums may be done explicitly, and we find 

F2a = |g| Yl ^(a;")^(x'")(tr([/"^)tr(t/")tr(t/-™) + {m^-m}+ 

n,m>0 

{n —n} + {m, n — m, —n}) 



(3.16) 

where we have left the divergent sums (which will not be relevant for the computation we 
are doing in this paper) explicit. The function F is related to the single particle partition 
function zy by 

/•OO 

F{e-^)= / dazvie-''), (3.17) 



Jp 

or explicitly 

2x 



F{x)=2\n{l-x) + - ^. (3.18) 

(1 — x)^ 

Thus, the contribution to Ci,i,_2 from F2a is 

C2a = ^(F(x)F(x) + 2F{x)F{x^)). (3.19) 

For the other two cases we could not compute the sums explicitly, but it is not difficult 
to numerically evaluate the desired coefficient at the transition temperature Xc = 2 — 
We find that the contributions of the three diagrams to Ci,i,_2 are given by 

C2a = 6.53536 X 10"^, 

C2b = -22.87088 x 10"*, (3.20) 
C2c = 9.16396 X 10"*, 

so that the total coefficient is 

Ci,i,_2 = -7.1716 X 10"* . (3.21) 

3.4- Three Loops 

The leading contribution to the coefficient Ci^i^_i^_i in (|2.8| ) comes from the fourteen 
three-loop diagrams of figure 2. We first give the expressions for each diagram in terms of 
propagators and spherical harmonic integrals. For diagram 3a, we find 

771 _ 9ym 

dt A^-^(0,aa6)Aj-^(t, aca)Aj,(t, Q!ac)(Aj^(0, aad) + Aj^(0, ttdc))- 

(3.22) 
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For diagram 3b, we find 

_ - _ - DiPPD^pp 

JpyJp + ^) 

{Aj^ (0, aab){DrAj^ (0, aad) + DrAj^ (0, Q!dc)) 

J dt (DrAj^ (t, Oac) Aj, (t, Oca) " Aj^ (t, aoc)^T Aj^ (t, ttca)) 

+ Aj„(0,Q!a6)(Aj-^(0,Q!ad) + Aj^{0,adc)) 

J dt {{ji3 + l)'A,^(t, a„c)A,,(t, Clea) " £'rA,^(t, a„c)£'r A,, (t, aea))| • 

(3.23) 

For diagram 3c, the result is 

" " 2 Ja(ja + 2)jp(jp + 2) 
y" dt(Aj^(0,aad) + Aj-^(0,adc)) 

+ + 1)' + (i7 + 1)'(J<5 + 1)')A,J0, aa6)A,,(t, aac)A,^(t, «ej 

- (j> + ifDrAj^it, aaa){4:DrAj^{0, aab)Aj,{t, aac) + 2A,J0, aab)DrAj,{t, aac)) 
-2{js + lfAj^{0,aab)Aj,{0,<^ac)} 

+ J dt{Dr Aj^{0, aad)+DrAj^{0,adc)) 

{ {j^ + ifAj^ {t, aca) {^DrAj, {t, aac)Aj^ (0, aab) + 2Aj, {t, aac)DrAj^ (0, aab)) 

- 2DrAj^ (0, aab)DrAj^ {t, aac)DrAj^ {t 

1 C^ca) 

-2{DrAj^ (0, aab)Aj^ (0, aac) + 2A,„ (0, aab)DrAj^ (0, aac))} ■ 

(3.24) 

For diagram 3d, we find 

_ _ 9yM jja-fX j-)f36X jja^p jjflSp 

J dtAj^ (t, aab)Aj^ {t, abc) {2Aj^ {t, Ocd) A^-^ {t, ado) + Aj^ {t, Ocd) Aj^ (t, a^a)) 

9yM_ j-^aj3X j-^-ySX j-^j^p j-^aSp 

4 

J dtAj^ {t, aab)Aj^ {t, abc) {Aj^ {t, acd)Aj^ {t, ada) - 4Aj-^ {t, acd)Aj^ {t, ada))- 

(3.25) 
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Evaluating diagram 3e, we obtain 

JpKJp + 2) 

/(3 26) 

Diagram 3f gives 

" ^ jx(jx + 2)jp(jp + 2) 

[4Aj^ (0, aab)Dr Aj, (0, a^d)DrA,^ (0, a^a) 

- 2((j^ + If + (j5 + 1)2) A,- JO, aab)A,,(0, acd)A,-jO, ad«) 

+ y" { Aj^ (t, Oab) (t, Ctftc) Aj^ (t, Ocd) (t, ada) 

(i, + l)2((J, + l)2 + (J^ + l)2) 

+ 2DrAj^ {t, aab)DrAj^ {t, abc)DrAj^ {t, acd)DrAj^ {t, aaa) 
-4(j^ + lfAj^{t, aab)Aj^{t, abc)DrAj^{t, acd)DrAj^{t, ada)}] 
_ 9ym D^'I^^D^^^D^^pD^'^p (3.27) 

[ADrAj^ (0, aab)DrAj^ (0, abc)Aj^ (0, a^a) 

- 2(j« + l)2A,jO,a„,)A,jO,MA,,(0, 

- 2A,-^ (0, aab)DrAj^ (0, abc)DrAj^ (0, 

+ y" {-Dr Aj„ (t, aab)DTAj^ {t, abc)DrAjg {t, acdjD^Aj^ {t, ada) 

+ Ua + l)^(j5 + l)^Aj^(t, aab)Aj^{t, abc)Aj^{t, acd)Aj^{t, ada) 
+ 2{js + 1)^ Aj^ (t, aab)DrAj^ {t, abc)Aj^ {t, acd)DrAj^ {t, ada) 
-4:{js + l)'^DrAj^ {t, aab)DrAj^ {t, abc)Aj^ {t, acd)Aj^ {t, ada)]] ■ 

For diagram 3g, we find 



F3g =gYME''^^E^^^{D^^^D^^^ - -D^^^D^^^ - -D^^^D^^^) 

J dtdt' Ajp{t' , ada)Aj^{t', acd)Ajg{t, abc)Aj^{t, aab)Aj^{t' - t, aac), 
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(3.28) 



where E was defined in the previous subsection. Diagram 3h evaluates to 

Fsh =9yM ■ D"7A^/3^Aga7pg^^p- 
jAUA + 2) 

j dtidt2 DrAj^{ti,aab)'^j^{ti,abc)Aj^{ti -t2,aca) 

(Ajg(t2, acd)DrAj^(t2, ado) - D-^Ajg{t2, acd)Aj^{t2, ado)) 

J dtidt2 DrAj^{ti,aab)Aj^{t2,ada)Aj^{ti -t2,aca) 

{Ajg {t2, acd)DrAj^ {ti, ate) - DrAj^ {t2, acd)Aj^ {h, abc))- 

Diagram 3i gives 

9yM pa.f3ppia.aj3 paS^ p5p^ 

j dtidt2dt3Aj^ {ti - t2, aab)Aj^ {ti - t2, oibc)Aj^^ (ts, acd) 

Ajsits, ada)Aj^{ti - ts, aca)Aj^{t2, aac)- 

For diagram 3j, we find 

Fsj ^g^M{D'^P^Df^P~^ - DPP^D^P~^)E^rap^af 

J dtdt'Aj^ (0, aab)Aj^ {t - t', aac)Aj^ {t', aad)Aj^ {t, aca)Aj^ {t', adc)- 

Diagram 3k yields 

Fsk =9ym ■ ^ D^P>^D''^~^&P-E'^--P 
Jx{jx + ^) 

J dtidt2 Aj^ {ti - t2, acd)Aj^ {ti - t2, ada) 

{2D, (0, aab)Aj^ {h, a„,)D, A,^ (^2 
+ Aj^ (0, aab)Dr Aj^ {h, aac)Dr Aj^ (ts 

-{ja + ^fAj^{0,aab)Aj^{ti,aac)Aj^{t2,aca)} ■ 

For diagram 31, we obtain 

/r„, = _ E°'^'^ R^7P E^ota ppaf 

J dtidt2dUAj^{t2 - ts, aab)Aj^{t3 - ti, aac)Aj^{ti - t2, aad) 
Ajp (ti, adc)Aj^ {t2, abd)Aj^ (ts, acb)- 
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For diagram 3m, we find 

(Aj^(t, adc)DrAjg{t, aad) - DrAj_^{t, adc)Aj^{t, aad)) 
Jx{jx + 2)Uj, + 2) 

+ 2L>^ A^-„ (t, aab)Aj^ {t, abd)Aj^ {t, ada)Dr Aj^ (0, a^a) 
-{js + lfAj^{t, aab)Aj^{t, abd)Aj^{t, ada)Aj^{0,aca)} ■ 
Finally, diagram 3n gives the result 

Tiajp T)B^a ( (~tpaXr~tXBa (~tp]3X (~tXaa 

F^^ = g^^^^—JL^ ^ 3 ^, , + 3 ^, , + D^'PDf'^ + D^'^D^^P 

"^^^ Jp Up + 2) + 2)\ jx [jx + 2) jx [jx + 2) 

{DrAj^ (0, a,b)Aj^ (0, aac)DrAj^ (0, aad) 

+ 2DrAj^ (0, a,b)DrA^^ (0, aac)A,^ (0, aad) 

-{j^ + 1)2A,„ (0, a,b)Aj^ (0, «ac) A,^ (0, } 

"^^"^JpOp + 2),.(,. + 2) I^S.0. + 2)+^ 
{L>,A,J0, a„b)A,^(0, a„,)L>,A,^(0, 

+ 2D, A,^ (0, aab)D^Aj^ (0, Oac) A,^ (0, a„d) 
+ 1)2A,„ (0, ttafc) A,^ (0, a^c) A,^ (0, aad) } 

-^^^^j,0, + 2)j.(,. + 2) l^',.0. + 2)+^ 
{2I),A,„ (0, aab)Dr Aj^ (0, aad) A,^ (0, a^^) 

+ 1)2 + (j^ + l)2)A,jO,a„6)A,^(0,aad)A,^(0,abe)} ■ 

(3.35) 

5.5. Evaluation of three- loop diagrams 

We now evaluate the coefficient Ci,i,_i,_i of /35r|.^|tr(C7)|^ in ( p.5|) arising from the 
three loop diagrams computed above. The expressions for the diagrams are of the general 
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form: 

F=Y1 S'{j,m,e)I'{j), (3.36) 

where 5" is the term containing aU spherical harmonic factors (C, D, E) and /' is the term 
involving the propagators and the integrals over t. S' depends on the j's, m's and e's of 
the spherical harmonics but /' only on the j's. Using the identities of appendix B we can 
perform the sum over the m's and e's analytically and find S{i) = ^ e)- Then, 

the diagram can be written as a sum over j's 

F = Y,S{j)I'{j). (3.37) 

j 

We then expand in powers of tr(C/"^) and determine /(j), the coefficient of |tr(C7)|^ in 
this expansion. The contribution to h from this diagram may then be written in the form 

Ci,i,_i,_i = ^^(j)/(j). (3.38) 

j 

We can find S for any diagram by using the identities given in appendix B, and the 
corresponding integral / can be found in the Mathematica file The sum ( p.38|) over j's 
is the only part of the calculation that we have generally had to perform numerically. 

Diagram 3a is the simplest and we can actually calculate it analytically. The expres- 
sion ( |3.22| ) above may be simplified using the angular momentum sums of appendix B. We 
find 



^3a ^ \_ V a(a + 2)6(6 + 2)c(c + 2) ^,+,+2 

^1,1,-1,-1 I87r4 ^ (a + l)(6+l)(c+1^3 



~c l^^-^n^^-^n^^-^J (3.39) 
{ (x"+^ + 1) (x"+^ + 2) - (c + 1) ln(x) (2x^^+2 + 3x"+^) } . 

Note that this sum contains a logarithmic divergence in the sum over c when taking the 
x-independent term (equal to 2) in the curly brackets. 

Such logarithmic divergences appear in the contribution to Ci^i^_i^_i from most of our 
3-loop diagrams, since there are terms where only two of the three unconstrained sums over 
loop momenta have exponential damping factors. These divergences arise from non-planar 
1-loop subdiagrams, as we explain in the next section when we discuss the corresponding 
counterterms. 



22 



In practice, when doing the computation we will place a cutoff J on the angular 
momentum sums. We can then write the result with a cutoff in the form 

= « x: ^ + + ^(v^), (3.40) 

c=l 

where a is the coefficient of the logarithmic divergence, which we can define by 

a = lim (J + l)(Ci,i,_i,_i(J) - J - 1)), (3.41) 

J — >oo 

and! 



c=l 

In particular, from the expression above for the 3a sum ( p.39|) we find 

«3a = ((13^ + ^°S(1 - = -4-0356 X 10-^ Ci^"-/-;^-^ = -3.666 X 10-^ 

(3.43) 

where we have evaluated ctsa and (^^a./j^^^e^ ^^^q transition point Xc = 2 — -\/3 0. 

For the 13 remaining diagrams, we were unable to evaluate the sum over j in ( |3.38| ) 
analytically. However, we have performed the sum numerically (at the phase transition 
temperature), and we will present the results in section 3.7. 

3.6. Counterterm diagrams 

We have seen that the contributions to Ci^i^_i^_i arising from individual 3- loop dia- 
grams contain logarithmic divergences. It turns out that all of these must cancel in the sum 
over diagrams. To see this, note first that no single-trace counterterm can contribute to 
the coefficient of |tr(t/)|*, since such a contribution requires four index loops, while planar 
counterterm diagrams at order are two loop diagrams with only three index loops. In 
fact, the only possible counterterm contribution comes from double-trace counterterms0 
of the form ti{AA)tY[AA) , which give rise to diagrams of the form shown in figure 3. But 
no such counterterm is gauge invariant, so in any gauge invariant regularization scheme 



^ To improve numerical convergence we use Aitken's method : if lim„^oo an = r then 
hmn^oo — ^, — o""*"^ = r, but the convergence of the second sequence is faster. 

See, for instance, [10,11] for recent discussions of how double-trace terms can contribute at 



leading order in the large A'^ limit. 
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(such as dimensional regularization) , there are no counterterm contributions at all, and 
all divergences must cancel in the sum over diagrams. Our regularization scheme does not 
respect gauge invariance, but the coefficients of logarithmic divergences are insensitive to 
the regularization scheme, so we must still find that all divergences cancel.0 

The preceding argument does not mean that we can ignore counterterms altogether. 
Indeed, the finite contribution resulting from the sum over diagrams does depend on the 
regularization scheme, and it is crucial to include the contributions of finite counterterms 
in order to obtain the correct result. The counterterms that can contribute take the form 

LcT = (^^^ ^ {citr{A,A,)tr{AjAj) + 2c2tr(A,A,)tr(A,A,)) , (3.44) 

where ci and C2 will be finite numerical coefficients chosen so that the results of calculations 
in our scheme match with results using dimensional regularization. 



Figure 3: Counterterm contribution to the free energy at order A^. 



At order A^, this vertex contributes to the free energy via counterterm graphs of the 
form depicted in figure 3. Following the same procedure as in the previous subsections, 
it is not difficult to isolate the coefficient Cc of (7y^j|tr(f/)^| in the contribution of this 
diagram to the free energy : 



,+2 



(3.45) 



Using spherical harmonic identities from appendix B and performing the sum we find: 

2 



2 

TC 



4 ^ Ci + -C2 



X 



{1-xy 



+ ln(l - x) 



which may be conveniently written as 



Cc = -(18ci + 12c2)a3c 



(3.46) 



(3.47) 



In comparing with dimensional regularization, the coefficient of 1/e poles will be proportional 
to the coefficient of the logarithmic divergence in a cutoff scheme. 
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In order to determine the counterterm coefficients ci and C2, it is enougli to look at 
the simplest correlator to which the counterterms ( p.44|) contribute, namely the nonplanar 
one-loop four point function on IR^ with two external legs attached to each of the index 
loops. As we describe in detail in section 4.4 and appendix A. 5 below, the coefficients 
are determined by demanding that the combination of logarithmically divergent one-loop 
diagrams (shown in figure 6) with the counterterm diagrams reproduces the result for the 
same correlator evaluated in dimensional regularization. For our regularization scheme 
with a sharp cutoff, we find 

ci = C2 = ^, (3.48) 

so that 

Cc = -^«3a. (3.49) 

This net contribution from the counterterm diagrams is all we need to complete our 
calculation. However, it turns out that a very useful check of our three loop results arises 
from splitting up the counterterm contribution into pieces associated with the individual 
one loop diagrams of figure 6. Note that each of these is logarithmically divergent, and 
appears as a nonplanar subdiagram in exactly one of the three loop vacuum diagrams 
of figure 2 (obtained by joining up the four free vector lines of the one loop diagram in 
pairs in such a way that we obtain a planar three-loop diagram, and replacing any a or 
c lines with the corresponding effective vertex). Thus, the logarithmic divergences in the 
contribution to h from three loop vacuum diagrams are directly related to logarithmic 
divergences in a specific set of non-planar one-loop four point diagrams.li^l If we denote by 
L^rp the counterterm Lagrangian density necessary to make this particular set of one-loop 
diagrams agree with dimensional regularization, then the combination of the three loop 
diagram X and the counterterm diagrams associated with L^j, should be finite, providing 
a check on the divergent part of each individual 3-loop diagram. 

In addition to the overall result ( p.49| ) we therefore define partial contributions 
arising from the counterterm diagrams associated with a given three loop diagram X. 



Tlie coefficient of |tr([/)|^ has no additional divergences from planar one-loop subdiagrams 
of tfie tfiree-loop diagrams, since we have seen that there are no single-trace counterterms that 
could cancel them (or remove the regulator dependence of the finite results if any such divergences 
cancelled upon summing diagrams). 
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These will generally have both a logarithmically divergent piece, and a finite piece, and 
take the forniEl 



C5 = «5lnf 1 +C. 



X finite 



(3.50) 



The values of etc and C^"^**^ for each diagram are tabulated in the next section together 
with our results for the three loop diagrams. 



3.7. Results 

In this section, we tabulate our numerical results for the contributions to Ci,i,_i,_i. 
For each three-loop diagram, we give the coefficient a of the logarithmic divergence and 
the remaining finite piece, as defined in ( |3.41| ) and (|3.42| ). We also give the coefficients etc 
and d™^^ for the logarithmically divergent and finite parts of the associated counterterm 
diagrams, as defined in (|3.5C1| ). We find: 



1 1 1 O ITT. O T>T 

uiagi am 


Oi/ CK3q, 


/^finite 

-1 -1 


Q^c/ CK3a 


^finite / 


3a 


1 


-3.666 X 10"^ 


-1 


-107/120 


3b 


-1 


-4.805 X 10"^ 


1 


13/60 


3c 


-3/4 


2.224 X 10-5 


3/4 


1/20 


3d 


5/4 


2.21 X 10-5 


-5/4 


-47/60 


3e 


1/2 


-5.85 X 10-5 


-1/2 


-7/30 


3f 


-9/4 


-2.69 X 10-^ 


9/4 


3/5 


3g 


-3 


1.15 X 10-^ 


3 


1/5 


3h 





1.91 X 10-^ 





-1/10 


3i 


5/2 


-6.6 X 10-^ 


-5/2 


1/12 


3j 


-3 


2.414 X 10-4 


3 


11/10 


3k 


1/2 


1.42 X 10-4 


-1/2 


7/60 


31 


5/4 


-1. X 10-4 


-5/4 


1/24 


3m 


2 


4.6 X 10-5 


-2 


1/15 


3n 


1 


-9.2 X 10-5 


-1 


-29/30 


Sum 





-4.5 X 10-4 





-1/2 



The constants .^1,^2,^3,^4 appearing in Appendix A all evaluate to the same value with 
the choice of damping function (essentially a step function) used here, while the UV cutoff M is 
simply J/R. 
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The numerical results in the middle column of the table are all accurate at least within 
the number of digits appearing in the table. We estimate the maximal total numerical error 
of our result for the sum to be less than 3 x 10~^. This accuracy could be improved with 
additional computer time, but we found no reason to do this since we are only interested in 
the sign of b. Note that the sum of all divergent contributions vanishes (without including 
counterterms) . 

Note that, in the table above, + = for every diagram X; that is, the sum 
of any three-loop vacuum graph and its associated counterterm graphs is finite. This sum 
depends on the renormalization scale fi (see ( |3.42| ) and ( |3.50| )) but should be independent of 
the regulating function R{q) 0. This /U-dependence vanishes upon summing all diagramsEl, 
leaving us with the scheme-independent answer 

Ci,i,_i,_i = -4.5 X 10"^ - «3a/2 = -4.3 X 10"^ (3.51) 

3.8. Effective potential 

We can now put together our results for the terms of interest in our effective potential 



( p78| ) for the eigenvalues, evaluated at the deconfinement temperature: 

^12 = 4.8112 X 10"\ 
/3cCi,i,-2 = -9.4447 X 10"^, (3.52) 
= -5.7 X 10-^ 

From ( |2.11j ) we find that the coefficient b in the effective potential ( |2.10| ) for ui is 



b = f3cCi,i,- 



1,-1 



/U2 

4 1 QK/I XX 1 n-6 (3.53) 



~ -5.7 X 10"* - 1.854 X 10 



-5.7 X 10"^ 



In order to verify this independence - and as a check on the logic of our regularization scheme 
and our numerics - we have recomputed a/asa, Cf"*!:i^_i, Oc/asa and Cc^^^^/a^a for diagrams 
3a, 3b, 3c, 3g and 3n with a different regulating function (we took R{q/M) to be a double step 
function, R{x) — 1 for x < 1, R{x) = ^ for 1 < x < 2 and R{x) = for x > 2). We evaluated these 
quantities analytically for diagram 3a and numerically for all the other diagrams. As expected, in 
every case the coefficient of the logarithmic divergence a/asa and ac/cxsa was unchanged in this 
new regulating scheme. Also, as expected, while finite parts of each diagram and its associated 
counterterm yielded diflFerent values in this new regulating scheme, the sum of every diagram with 
its associated counterterm graphs was unmodified. 

This follows because T^r ,. = Y^r ,. a'^ = 0. 
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Note that b is the sum of two terms, one of which is manifestly negative. In ( |3.53|) , 
however, the dominant contribution came from the first term which, in principle, could 
have been either positive or negative. This suggests that there could exist also theories for 
which 6 > 0; we will discuss this further in the final section. 

Since we have, after a laborious calculation, determined that 6 < 0, we may conclude 
that the large deconfinement transition of pure Yang-Mills theory on a small is of 
first order. 

4. A Gauge Invariant Regularization Scheme on 

In the previous section we have computed a set of two and three loop Feynman dia- 
grams to determine a particular term [b in equation (|1.4|) ) in the Wilsonian effective action 
for finite temperature Yang Mills theory on S^. We found that while b is finite, individual 
diagrams that contribute to b diverge logarithmically. In order to obtain the finite physical 
value of b we needed to sum contributions from the various diagrams, at which point the 
logarithmic divergence cancels and we are left with the finite result of interest. The process 
of isolating a finite piece from the difference of divergent sub pieces is delicate, and will 
yield the correct answer only if the regularization procedure respects gauge invariance. In 
this section we will expand on the discussion of section 2.4 to describe in more detail the 
regularization procedure that we employ in our computation. 

As described above, computations of Yang Mills theory on are most simply per- 
formed in the Coulomb gauge diA"^ = 0, where the index i runs over the three spatial 
indices of the 5"^. Further, we found it most convenient to regularize all diagrams by 
truncating the spherical harmonic sums at spherical harmonic number n (in fiat space this 
corresponds to imposing a hard momentum cut off at momentum E{n)/R where E{n) ~ n 
is the energy of the n*'^ spherical harmonic mode). This regularization scheme is not gauge 
invariant, but should yield gauge invariant results when employed with a bare Yang Mills 
action that includes an appropriate set of non-gauge-invariant counterterms. The appro- 
priate counterterms may, in principle, be uniquely determined (up to the usual ambiguity 
in the definition of the Yang Mills coupling constant) by demanding that correlation func- 
tions computed by this theory obey the Ward identities that follow from gauge invariance, 
together with local Lorentz invariance. 

In this section we will explicitly determine some of the counterterms that will render 
our non-gauge-invariant regularization scheme effectively gauge-invariant. These counter- 
terms fall into two classes; counterterms that would be needed even in flat space, and 
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counterterms that are proportional to the spacetime curvature. It wiU turn out that no 
counterterm of the second type (those proportional to spacetime curvature) contributes 
to the computation of b, so we will content ourselves with determining only those coun- 
terterms that appear even in flat space. These counterterms may be determined rather 
simply by choosing them to ensure that certain Green's functions (following 't Hooft we 
use the A^{p)At,{—p) two-point function, as well as a four-point function of gauge fields) 
agree with the same Green's functions evaluated using dimensional regularizationlll. 

In §4.1 we explain our regularization method in detail and give a simple example of 
how it works. Sections 4.2 and 4.3 are devoted to tests of the validity of our regularization 
scheme. In §4.2 we determine all quadratic counterterms of the first type (those that appear 
in flat space) to order A. Even though these counterterms do not actually contribute to 
our main computation in this paper, we use them to test the validity of our regularization 
procedure. First, we verify in §4.3 that our results are consistent with the Slavnov- Taylor 
identity and with Lorentz invariance. As another test, in appendix A. 3 we verify that our 
results lead to the correct free energy at infinite volume to order A. Finally, in §4.4 we 
proceed to use the same methods to compute the counterterms that we actually need for 
the computation of b; these are a set of double-trace counterterms at order A^. 

4.1. General discussion and a simple example 

The regularization we will analyze in this section is a slightly more general regular- 
ization scheme than the sharp cutoff which was used in the computation of the previous 
section. We include damping functions R{>/q^ /M) and R{qo/A) for the momentum of 
each internal Ai line of a given diagram, and damping functions R{^/(p / M) and R{qo/A) 
{q^ = q^qi, i = 1,2,3, and we take A>A, M>M, A>M and A > M for conve- 
nience) for the momentum of each internal Aq or ghost line. These functions are chosen 
so that i?(0) = i?(0) = ^(0) = ^(0) = 1, i?'(0) = i?'(0) = ^'(0) = ^'(0) = 0, and 



More precisely, we compare with a form of dimensional regularization that is tailored to deal 
with Yang Mills theory in the Coulomb gauge. This so-called split dimensional regularization 



scheme |12| , |13| , |14] , |15| separately extends the number of dimensions participating in the Coulomb 
gauge condition (from 3 to 3 — e) and the number of other dimensions (from 1 to 1 — e'). One may 
worry that this regularization procedure is not Lorentz invariant (since integrals over temporal and 
spatial momenta end up being regulated differently), however the breaking of Lorentz invariance 
really comes from our choice of gauge rather than the regularization scheme. In practice, we will 
apply split dimensional regularization in the limit where e' 0. 
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R{x oo) = R{x oo) = R{x oo) = R{x oo) = 0. We choose to treat Aq and 
ghost hnes differently from Ai hnes because, in our calculation of 6, it was convenient to 
integrate out Aq and the ghosts directly in the action. This can be done with no regu- 
larization subtleties in diagrams with both Ao/ghost lines and Ai lines, provided that we 
take the scale M ^ M and A ^ A. We will sometimes be lazy with our notation and 
write R{^/M) {R{^ / M)) as R{q/M) {R{q/M)). 

In perturbation theory, correlation functions are obtained by evaluating all contribut- 
ing Feynman diagrams, each of which may be written as an integral over internal momenta. 
In general these integrals diverge and must be regulated; in this section we will explain 
how one may convert the simple minded regularization scheme described in the previous 
paragraph into dimensional regularization by an appropriate choice of non gauge invariant 
counterterms. In the rest of this subsection we will demonstrate our method on a 'toy' 
regularized integral0 

d^qdqo Ri^)R{^) 

(2vr)4 q2 + qi ' ^^''^ 
and its counterpart in split dimensional regularization (SDR) 

d^qdqo 1 
SDR (27r)4 q^ + ql 



(4.2) 



The go integral in ( |4.1[ ), ( |4.2| ) is finite (which is why we ignored the R regulators in ( |4.1[ )) 
and may easily be done to yield (from now on we suppress explicit reference to the regulator 
in intermediate steps) 

iDTr-^ J q 



In split dimensional regularization (|4.3|) evaluates to zero, whereas in the damping function 



regularization scheme (|4.1|) evaluates 



2 2 

a{p) = M^C2 + - (4.4) 



This integral is slightly different in form from those that will appear in our actual expressions 
below; in particular the integrand contains a single propagator but two copies of the regulator 
function. Thus, it should merely be thought of as a simple divergent integral that illustrates all 
the complications that arise in the actual process of regularization. 

The finite piece arises because JJ^ dqR{q)R'{q) = —\- 
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where we have defined 



Co 



1 



47r2 
1 

4^ 



dqqR{qf, 
dqqR{q) R"{q). 



(4.5) 



Thus, in order that our damping scheme agrees with dimensional regularization, we should 
introduce a counterterm that contributes —a{p). 

Not every integral we encounter will be simple enough to explicitly evaluate in split 
dimensional regularization (as / of (O) was). However, it will always be possible to de- 
compose the integrals of interest into the sum of a complicated but convergent integral and 
an easily evaluated divergent integral; this will be sufficient to determine the correspond- 
ing counterterms. As an illustration, we reevaluate the counterterm corresponding to the 
integral / in a perversely convoluted manner. Under the change of variables q p~q (an 
allowed variable change in both regularization schemes), / becomes 

d^q 1 



(27r)4 (p-g)2 + (p„_g,)2' 

Performing the go integration as above yields 



1 



167r3 



d^q 



1 



\/{p-(iy 



167r3 



q-p p^ 



q^ 



2r 



+ 



2g4 



+ 



(4.6) 



(4.7) 



We may now rewrite (|4.6| ) as a sum over a manifestly convergent piece and an easily 
evaluated divergent piece a^ 



1 



167r3 
1 

16^ 



d^q 



\/{p- 



d^q 



\/{p-(lf 



1 + 



+ 



1 



167r3 



d^q 



1 + 



q-p 



p^ 



q-p 



+ 



p^ 



2(g2 + a2; 
?,{q-pf 



+ 



"iiq-pf 



2q^{q^ + 



(4.^ 



2(g2 + a2) 2q^{q^ + a- 



In (split) dimensional regularization with d = ?> — e, the second line evaluates to 



327r3 



d^q 



q^{q^ + a2) 



?>{q-pf 



P 



p^ 



lim 

.2 



d'^q 



= lim 

P 



p- 



327r3 

2 



/ 4n 



\ e 



(4.9) 



247r2' 



19 



We introduce a parameter a in order to avoid artificially introducing IR singularities. 
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On the other hand, in the cut off scheme it evaluates to 

M'Ca + ^F2, (4.10) 

and ( [4.101 ) is equal to ( [4.9|) upon the addition of the counterterm —a{p). 

In this particular example, the convergent part of ([4.8| ) is easy to evaluate (it is equal 

_ 2 

to 24^)' t)ut this is not true in more complicated examples, and was not needed in order 
to evaluate the counterterm. 

Of course, the separation of (|4.8|) into a divergent and a convergent part is ambiguous. 
In our work ahead we will find it convenient to fix this ambiguity by demanding that the 
divergent piece (which we call the regulator dependent piece below) should evaluate to zero 
in split dimensional regularizationS. With this convention, the counterterm associated 
with any diagram is simply minus the regulator dependent integral evaluated in the cutoff 
regulator scheme. 

4-2. Single-trace quadratic flat- space counterterms at order \ 

As an explicit example, we can now proceed to compute the regulator-dependent piece 
of the gauge boson self-energy II^i, = —^{A^Aj^) . To set our conventions, we write the 
Yang-Mills action as 

S = \J d^xtriF^^F^^"). (4.11) 
The momentum space Coulomb gauge {diAi = 0) propagators take the form 

{AfAf) = 5''''5'- ( tEl2^^^\ , («) = ( \\ . (4.12) 

In addition, the gauge fixing procedure introduces a set of complex adjoint ghosts c, c with 
Lagrangian 

Cghost = -tr(ca^Ac), (4.13) 

where Di = di — igYM[Ai,*] is the gauge covariant derivative. The propagator for the 
ghosts is identical to that of the Aq fields. 



With this convention, the regulator dependent piece in the example of the previous para- 
graphs is the second term in (|4.8| ) minus 
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In the computations of this section we do not exphcitly integrate out Aq and c as we did 
in the previous section; of course this does not affect the results. 

In Appendix A.l we define several regulator-dependent constants and functions of 
external momentum that arise in our calculation. In Appendix A. 2 we depict the dia- 
grams contributing to the gauge boson self-energy, and list our results for their regulator- 
dependent contributions (defined above). Adding these together yields the following result 
for the regulator-dependent contribution to the gauge boson self-energy (see Appendix A.l 
for notation) : 

J^ik^^ ^ 8^ ("^) ~ (p'^ +PQ)9ik] - ^ {p'^gik + ^PiPk) F2 

+ iTZJiPomk In -4 -f -—rp g,k\n\-^\+ ^:r^PrPk In 



247^2^^^-"^ V ^4 y ' 407r2^^*'^ V^l; 1207r2^*^'^^" V^l' 

\^r--^P.P.^{^).^P.^^ (.1. 

in^r^ = --iL.^ ln(^]- ^ R - J^.^ + 2M^ (C. - C.) + -J-..^ In ( J 



_ In _ - + 2M' (Ci - C2) + :tp' In , 

(4.17) 

From the discussion above, the required quadratic counterterm Lagrangian must be 
chosen to precisely cancel these regulator dependent contributions, in order to give agree- 
ment with dimensional regularization. Thus, we must have 

£et = -tr(A'^A-)n(,^^). (4.18) 

In the next subsection, we perform two consistency checks on these results. 

4.3. The Slavnov- Taylor identity and SO {A) invariance 

In this subsection, we first use the fact that our result for the self-energy must be 
consistent with gauge invariance and 5*0(4) symmetry at short distances to demonstrate 
the consistency of our results for the logarithmic divergences in ( ^.15| )- ( [4.17| ). In particular. 
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since the logarithmically divergent and finite contributions must satisfy various Slavnov- 
Taylor identities for these symmetries independently, we can determine the structure of 
the former without any knowledge of the latter. 

We begin by considering the Slavnov- Taylor identity relevant for the gauge symmetry 
in Coulomb gauge. As usual, we start with the Euclidean gauge-fixed action 

S = Jd^x !^Cym + Y^i^iA'f - cViAcj (4.19) 

and take e ^ to get the Coulomb gauge. The BRST charge Q satisfies 

[g, A^] = D^c, 

{Q,c} = ic'^, (420) 
{g,c} = ^ViAi. 

To obtain the Slavnov- Taylor identity, we study the partition function with sources added 
for the operators generated by the BRST transformation 

Z[J^,C,K^,L] = J exp|-5 + Jd'^x {J^A^ + Cc + cC + K^[Q,A^]-L{Q,c})Y 

(4.21) 

Performing the change of variables A ^ A + [eQ, A], c ^ c + [eQ, c], c ^ c + [eQ, c] in the 
path integral, we eventually obtain the standard identity 

+ T777 = 0' (4-22) 



6Ki^ 5c 5L 



where F is the IPI efi"ective action less the gauge fixing term. From this, we may easily 
derive a 1-loop Slavnov-Taylor identity relating the self-energy, H/^u, and the coefficient 
$^ of the K^c term of F: 

d^U^, + {-d^g^u + d^d,)^i, = 0. (4.23) 

An analogous relation arising from SO {4) invariance is difficult to obtain since we are 
working in a noncovariant gauge. Fortunately, a simple restriction that arises by requiring 
50 (4) -invariance of the S- matrix will be enough for our purposes. The specific condition 
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that we will impose is the existence of a double pole in the full propagator at zero momen- 
tum. This requirement, combined with the weaker Slavnov- Taylor identity d^d^M^j^ = 0, 
restricts the local part of Il^jy to be of the form 

liij =C {[p'^ + pl]gij - piPj) 

n^o = Dp.pQ (4.24) 
noo = -(C + 2I^)p2 

where C and D are dimensionless constants. It is easy to demonstrate, using the definition 
of that the full Slavnov- Taylor identity ( ^4.23| ) fixes C + D, leaving one degree of 



freedom that is in principle determined by a wave function renormalization condition. The 
logarithmic divergences must have this structure independent of the finite contribution to 
the local part of li^y. That our result ( [4.15| )- ( [4.17| ) is consistent with these conditions is 
easy to verify. 

As a second, and less formal consistency check, we have used our regularization scheme, 
together with the counterterms (|4.18|) , to compute a physical quantity; the two-loop free 



energy of Yang-Mills theory at infinite volumellll. The counterterms computed in the 
previous subsection play a crucial role in our calculation, which we present in detail in 
appendix A. 3. Our final answer, F2-100P = VXT'^/72, agrees with the previously computed 
result (using dimensional regularization in Feynman gauge) . We regard this as a rather 
nontrivial check of our regularization scheme. 

4.4- The tr{AiAj)tr(AkAi) counterterms 

We have seen in section 3.6 that the counterterms required to evaluate b take the form 
of double-trace terms quartic in the spatial components of the gauge field. We will now 
follow the method of the previous subsection (requiring the order }? /N"^ contribution to 
the four point function {AiAjAkAi) to agree with the result obtained by (split) dimensional 
regularization) to evaluate the coefficients of the two possible counterterms of this form, 
given in ( |3.44|) . 

The one loop diagrams contributing to the nonplanar part of the four-point correlator 
are depicted in figure 6 in appendix A. 5. It follows from power counting that the leading 
divergence in each of these diagrams is logarithmic. As a consequence, the regulator 



As far as we are aware, this is the first time that this computation has been done in Coulomb 
gauge using any regulating scheme. 
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dependent part of each of these diagrams may be evaluated with aU external momenta set 
to zero and has a divergent part proportional to the single integral 

(4.25) 



In Appendix A. 5 we list the coefficient of this divergent integral computed for each of the 
diagrams with a particular index structureii. To obtain the full expression, we must also 
sum over distinct permutations of indices. We list the results -R^*^; (for the coefficient of 
the integral ( [4.25|) ) , as well as the corresponding contribution to the counterterm, diagram 
by diagram in Appendix A. 5. 

Summing over the expressions in Appendix A. 5 we find that the sum of the diagrams 
depicted in Appendix A evaluates to 



1_ 



d{d + 2) 



Notice that the first term in ( [4.26| ) is simply zero in the cut off regulator scheme (on setting 
(i = 3). However, in dimensional regularization this term evaluates to 

+0(e^)l x- = -T#^- (4-27) 



It follows that perturbative computations in the damping function scheme must be accom- 
panied by the counterterm 

= (^) (125^ ('^^^) + 6^) [tr(^.^.)tr(A,A^) + 2tr(A.A,)tr(A.A,)] , 

(4.28) 

where the ^s are regulator dependent constants defined in appendix A.l. For the sharp 
cutoff used in section 3, all ^'s give the same result, so the term involving ^'s evaluates 
to zero. 



We keep the number of spatial dimensions, d, explicit (as opposed to setting d = 3), in order 
to determine the finite counterterm needed to bring our result into agreement with that of (split) 
dimensional regularization. 
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5. Validity of Perturbation Theory 

In this section we will determine the precise regime of validity of perturbation the- 
ory for pure Yang Mills theory on a sphere of radius R, in order to make sure that the 
computation we described above is valid. Naively, perturbation theory is good whenever 
^qcdR ^ I5 since the running coupling constant is then small at all scales above the 
scale 1/R oi the classical mass gap. In thermal Yang Mills theory it turns out that this 
expectation is modified by IR divergences; as we explain below, perturbation theory is 
valid at small AqcdR only for TR <^ a(T) ' where A(T) is the running 't Hooft coupling 
at the energy scale T. It follows in particular that for AqcdR ^ 1 perturbation theory is 
good at TR ~ 1, which is the regime of interest for this paper. 

5.1. Review of IR divergences in flat space 

In this subsection we review the well known effects of IR divergences on thermal Yang 



Mills theory in fiat space; see ||T^ and references therein for more details. 

Perturbative computations in Yang Mills theory on IR'^ x (where the 5"^ is a thermal 
circle) are beset by IR divergences, as is easily seen from power counting. IR divergences 
arise from the w = sector (w is the Euclidean energy) of the theory. Consider a Feynman 
diagram made up entirely of w = modes. Let q be the scale of spatial momenta in such 
a diagram. Each additional loop is accompanied by a factor of ^ (from vertices and 
propagators) times q'^ (from phase space), giving a net factor of Consequently, higher 
loop graphs are increasingly infrared divergent. 

These infrared divergences are cured by the fact that the gauge field A is effectively 
massive. Working in Feynman gauge, the one loop self energy of the Aq field at zero energy 
and momentum, noo(0, 0), is nonzero and of order AT^. As a consequence, Aq is effectively 



massive with mass of order niei = y noo(0, 0) ~ y/XT. Consequently, infrared divergences 
in loops involving Aq are cut off at this mass; thus, the effective loop counting parameter 
for Aq loops with = is ~ -\/A- So, Aq IR divergences change the perturbative 
expansion parameter from A to a/A. The first fractional power of A that appears in the 

3 

expansion of the free energy is A 2 (from a one loop graph using a mass corrected propagator 
for Aq). The next fractional power, A2, follows from the electric mass regularization of a 
3-loop IR divergence. 

IR divergences involving the spatial gauge field Ai are more serious. It turns out that 
the Ai self energy at zero momentum vanishes at one loop, but is nonvanishing at two 
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loops. As a consequence, the effective mass for Ai is of order rrimag ~ V A^T^ = TX. IR 
divergences involving spatial Ai fields are cut off by this mass; as a consequence the effective 
loop counting parameter for loops involving spatial Ai is of order one. Graphs of low 
enough order do not suffer from spatial IR divergences; however, a detailed investigation 
reveals that an infinite number of graphs contribute to the free energy at order A"^ and 
higher. In summary, the free energy may be expanded up to order A^; all coefficients in 
the expansion of the free energy to this order are perturbatively computable, and have 
been computed (see []T^ and references therein). Higher order terms are, in principle, 
inaccessible to perturbative analysis. 

The generation of an electric mass simply refiects the fact that the high temperature 
dynamics of Yang Mills theory deconfines. Indeed, space is filled with a plasma of charged 
particles of density ~ T'^. As each of these particles carries a charge V^, the screening 
length of this plasma is 1/ VXT, explaining the magnitude of rriei described above. 

The generation of a magnetic mass may be explained from the observation that Yang 
Mills theory on IR^ x reduces, at high temperatures, to a (Euclidean) 3 dimensional 
Yang Mills theory with an effective Yang Mills coupling constant AT, coupled to an adjoint 
scalar field of much larger mass Tat. The low energy dynamics of this theory is simply 
that of pure 3 dimensional Yang Mills theory, which non-perturbatively develops a mass 
gap of order AT. 

5.2. IR behaviour on 

We now turn to a study of the IR behaviour of Yang Mills theory on S^. Yang 
Mills theory on has a mass gap 1/R even classically. As a consequence, even ignoring 
the dynamical mass generation, the power counting arguments of the previous subsection 
indicate that (assuming TR ^ 1) loops of low energy Aq and Ai fields are both weighted 
by the effective coupling 

X^ff XTR^m^agR. (5.1) 

Perturbation theory is valid when this effective coupling is small. When rrimagR ~ 1 this 
effective coupling is of order one, and perturbation theory breaks down 11. 

In summary, finite temperature Yang Mills perturbation theory on is useful pro- 
vided both that AqcdR ^ 1 and that temperatures are low enough so that X{T)TR <^ 1. 



Once naive perturbation theory breaks down, the effects of dynamical mass generation are 
important. In this regime the correct way to proceed is to mimic the flat space analysis, and to 
shift the bare quadratic action by Afj_Ilf^^{0,0)A^, where the first index of IT refers to the energy 
and the second to the spherical harmonic number on S^. At low enough temperatures this eflFective 
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6. Conclusions 



In this paper we liave computed tiie leading perturbative correction to the thermal 
partition function of pure SU{N) Yang-Mills theory on around the phase transition 
point, and found that it leads to a first order deconfinement phase transition. The analy- 
sis, requiring diagrams with up to 3 loops, is quite complicated, and the interesting result is 
a single number ( p.53|) which governs the order of the phase transition. As described above, 
we have subjected our formalism to various consistency checks, including the cancellation 
of all divergent contributions, but we do not have any way to independently verify the cor- 
rectness of our final result ( |3.53| ) . It would be useful to have an independent computation 
of (|3.53| ) as a check of our results. 

In the pure Yang-Mills theory on we found that b is negative; a similar result 



was found in the corresponding analysis of various quantum mechanical systems [ll9|,P0 
It would be interesting to compute the sign of b in other 3 -f 1 dimensional (or lower 
dimensional) field theories, and to see how it depends on the matter content and on the 
various coupling constants of the theory. In particular, we are planning to compute the 
value of b in the 3 -f 1 dimensional Af = 4 supersymmetric Yang-Mills theory on S^, to see 
if the order of the deconfinement phase transition at weak coupling is the same as the first 



order behaviour found at strong coupling [21,22 



Since b is negative in all examples that have been analyzed so far, one might conjecture 
that for some reason it always has to be negative. However, it is easy to see that this is 
not the case, at least when one adds additional scalar fields with arbitrary couplings. For 
example, let us consider aO + 1 or 1-1-1 dimensional gauge theory with a massive adjoint 
scalar field $, with a potential of the form 

= tr(m2$2 + C4^7^M$' + cegt^M^' + C8^7^m$'), (6-1) 

where the Ci are kept fixed in the 't Hooft large limit. By analyzing the vacuum diagrams 
up to 3-loop order, it is easy to see that cq contributes linearly to b at order (through a 



mass is of order 1/R, while at higher temperatures the effective mass crosses over to its flat space 
value. For instance, the effective Ai mass is of order 1/R for XTR <C 1, but is given by rrimag for 
XTR 3> 1. Consequently, for T ^ 1/XR, the flat space analysis of the previous subsection applies, 
and perturbation theory breaks down. Were we to ignore the dynamically generated contribution 
to the masses, we would be faced with a paradox. The free energy would receive contributions 
proportional to increasingly high powers of R, in conflict with extensivity. 
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diagram similar to 3n of figure 2), while cg only contributes to b at higher orders. Thus, for 
given values of and C4, we can achieve any sign for the leading perturbative contribution 
to b just by varying cg. We may need to choose cq to be negative for this, but we can 
always choose cg to be large enough so that $ = is still the unique minimum of ( |6.1|) . 
In higher dimensions we have not yet been able to find a similar example involving purely 
single-trace interactions, but a potential term of the form c(7y^,^tr($^)^/A^ may be shown 
by similar arguments to lead to arbitrary values for b (already at order A). Thus, large 
weakly coupled deconfinement phase transitions may generally be of either first order or 
second order. 
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Appendix A. Details Related to Regular izat ion 

In this appendix, we provide many of the details behind calculations discussed in 
sections 3 and 4. 

A.l. Definitions of useful regulator- dependent constants and functions 

We start by presenting definitions for the various constants which encode the depen- 
dence of the results in section 4 on the regularization functions R, R, R and R. We also 
define three functions of external momentum that will arise in the expressions for individual 
diagrams presented later in this appendix. 



A* y Jo g2 + ^ 



M2 



27T 



oo 



^1 = 77- / dQo R{qo) 



00 
00 



1 f°° 
Bi = — dqo R{qo) 

I 

B2 = — dqo R{qo) 
27r J_oo 



00 
00 



Ci = 4^/ dqqR{q) (A.2) 







Cl = ^l^ dqqR{q) 

C^2 = ^y^ dqqRiqf 
1 f°° 

P^ = ^l dqqR{q)R"{q) 



Note the parameter jj, appearing in the definition of the constants Ai. This parameter 
is needed when splitting logarithmically divergent integrals in order to avoid introducing 
artificial IR divergences into the "regularized" pieces. Moreover, ji is identified with the 
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scale associated to (split) dimensional regularization when comparing results obtained in 
that scheme with those obtained in ours. 



A. 2. Diagram by diagram contribution to the regulator dependent piece of II 




SE3a SE3b 

Figure 4: Diagrams contributing to the gauge boson self-energy at 1-loop. 
Solid lines denote Ai propagators, dashed lines denote Aq propagators, and 
arrowed lines denote ghost propagators. 

We now list the regulator dependent contribution, III'^^'* , of each diagram in figure 4 
to the self-energy Yl^i, = —^{A^Ai,) at momentum p : 

^ fj,iRD)\ 2 2^ 2 9 2 31 \ fA2M\ 

1/2 . ^ r. 269 17 2 1 2 

-j-^{p 9ij + ^PiPj) F, + ^^P^PJ - ^P 9ij + Y^Po9^J, 

(A.4) 
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(A.5) 

X (nS'JO,,. = (p.[^2], - 2M(^2).,) , (A.6) 
\ (4S) , , = (M[^2]ii - P4^2]i) , (A.7) 



A 

A 



' {^%11)^=Im'C,9,j, (A.8) 



X (nS?/)^^. = 2AMB,Cmj, (A.9) 



£^))^^ = 2M^Ci. (A.14) 



v4.5. v4 check : the 2-loop free energy 



We now proceed to compute a physical quantity, the free energy of Yang Mills theory 
at order A, using our regularization scheme. We will find a result that agrees with the 
standard result obtained by traditional methods (utilizing dimensional regularization in 
Feynman gauge) . We view this agreement as a significant check on the consistency of our 
regularization scheme. 

The two-loop free energy of Yang Mills theory, at order A, receives contributions from 
six graphs depicted in figure 5. A set of one-loop counterterm graphs also contribute to 
the same order. We will find it useful to group each two-loop graph with a set of one- loop 
counterterm graphs, in the manner (and for the reasons) that we now explain. 

Consider the counterterm contributions to the free energy at order A. For the purpose 
of this discussion, it will be useful to regard counterterms that cancel different one-loop 
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FEld FEle FElf 

Figure 5: Diagrams contributing to the 2-loop free energy. 

contributions to the self energy (see above) as distinctil. Any one-loop contribution to 
the free energy with a counterterm insertion may uniquely be associated with one of the 
two-loop free energy graphs by "blowing up" the counterterm into its parent self-energy 
diagram. It is plausible (and true) that the sum of any particular 2-loop free energy 
graph with all its associated 1-loop counterterm graphs, is finiteS and independent of the 
damping function R{j^)- This fact suggests a natural grouping of graphs that we will use. 

Next, we list the results related to the diagrams, depicted in figure 5, that contribute 
to the 2-loop free energy. For each such diagram, we list the result of evaluating it in our 
damping function scheme, F^, the contribution arising from its associated one-loop self- 
energy diagram (see section 4 and the previous subsection), , and the result obtained 

by direct evaluation in (split) dimensional regularization, f!^^^^ (all divided by XV, where 

Furthermore, it is easy to convince oneself that only the counterterm with external vector 

lines yields a temperature-dependent contribution to the free energy. 

It may also seem natural to guess that this sum will equal the corresponding diagram evalu- 
ated in split dimensional regularization, but this is not precisely the case, as we explain in detail 
in the next subsection. 
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V is the volume of space) 
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aJ 225' 
dp 



2761 , 
216000' 



12877470 qHei/T-l)Jo p2(ep/T_i) 
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(A.16) 
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Ffeic = -2AMB2TJ2 
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Fi,%l = 2AMB,TY, 
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Po 



FEld 
FEld 



PO 



(27r)3 p'^{p'^+pI) 



d'p (pV^ -pV)(^2)^,(p) 

(27r)3 pHp''+pI) 

d^P [p^g'' -fp'){H2Mp) 



(A.17) 



(A.18) 
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piDR) 
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(A.19) 



108' 



^{CT) 
FElf 

ADR) 



-AMCiBiT^ 
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(A.20) 







^ FElf 

Summing the contributions of the individual diagrams, we find that the 2-loop free 
energy is 

(A.21) 



2 — loop 



72 ' 



in agreement with results previously computed using dimensional regularization in Feyn- 



man gauge 



2^ The result often seen in the literature is jj^T^, which differs from this by a factor of 2. 
The reason for this apparent difference arises from a choice of convention. Our definition of 
the propagators in (4.12) corresponds to a normalization of the basis matrices [t"^)'^^ such that 
{t'^)'^''{t"^)'"^ = S'^'^S''", which is equivalent to taking the quadratic Casimir of the fundamental 
representation as C{N) = 1. The convention most prevalent in the literature is C{N) = ^. This 
difference amounts to an effective difference in the definition of A which, when properly accounted 
for, gives an extra factor of 2. 

2^ We also find the same result for the free energy evaluated directly in Coulomb gauge using 
split dimensional regularization. 



46 



We may explicitly verify the diagram by diagram cancellation of divergences discussed 
above. In addition, we note that adding the contributions from a free energy diagram in 
our damping function scheme and its associated self-energy diagram does not lead to full 
diagram by diagram agreement with the result of (split) dimensional regularization. We 
discuss this in greater detail in the next subsection. 

A. 4- Diagram by diagram comparison with split dimensional regularization 

In this subsection, we will explain in detail how the sum of any particular two loop 
self energy graph plus all its associated 1-loop counterterm graphs differs from the same 
graph evaluated in dimensional regularization. 

If we consider summing a free energy diagram computed in our scheme with its cor- 
responding counterterm, this yields a result equivalent to evaluating the integral over the 
internal loop momentum from which the divergence arises before contracting the legs of 
this loop with the Ai propagator. After the evaluation, we perform this contraction but, 
at that point, the number of spatial dimensions is fixed at d = 3. On the other hand, when 
we compute the free energy diagram in pure (split) dimensional regularization, contrac- 
tion of the divergent momentum integral with the final propagator is done at unspecified 
d, leading to additional dependence on e = 3 — d and changing the finite result. 

As a check on our calculations, we now compute the difference between the finite parts 
of diagrams evaluated with our scheme and with (split) dimensional regularization. Any 
factors of d that arise only affect the evaluation of the logarithmic divergences and thus 
wc restrict attention to these. We write a generic logarithmically divergent integral that 
arises in the asymptotic expansion of a given self-energy diagram as 

[/(/, e)g^, + g{p\ e)p,p,] [ f^^ \ -. (A.22) 

J (27r)"^ ^qf2(q,2 _|_ 

Contracting this with an A^ propagator at momentum p, and evaluating the integral in 
dimensional regularization, we obtain 



l-ln(« 



(A.23) 



We see that the new finite contribution that arises due to the e in the (2 — e) factor is 
precisely — | times the coefficient of In// in the final result. This implies that, to get the 
(split) dimensional regularization result for diagram FEla (FElb), we should add ■^T'^ 
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( 




). It is easy to check that this is consistent with equations ( |A.15| ) and ( |A.16|) . 



Moreover, we see why this diagram by diagram discrepancy didn't ruin agreement in the 
final result for the free energy; the difference between our scheme and (split) dimensional 
regularization is proportional to the logarithmic divergence of a given diagram and all 
logarithmic divergences cancel when the diagrams are summed. 

A. 5. Coefficients of divergences in non-planar one-loop 4-point graphs 

In this subsection, we consider details related to the computation of the tr{AiAj)tr{AkAi) 
counterterm that is relevant for the computation of b in section 3. The relevant diagrams 
are depicted in figure 6. The seemingly chaotic ordering and labeling of diagrams was 
chosen for easy identification with 3-loop free energy diagrams via the correspondence 
discussed in the main text. 

We compute the counter-term by computing four-point functions of gauge fields, with 
the color index structure {tr{AiAj)tr{AkAi)) which arises from one-loop non-planar dia- 
grams. These diagrams generally have a logarithmic divergence proportional to 



We begin by listing the coefficient, of this logarithmically divergent integral for 

each diagram in figure 6 with the indices fixed as in the figure, computed in dimensional 
regularization with d — 3 — e : 




(A.24) 



M. 



(3a) 



[(5 + 2d)gabgcd + (5(i I9)gacgbd + (17 - 4:d'^)gadgbc\ + [d 



3)gabgcd 




abed 




(A.25) 

[((i^ - 2)gabgcd + {d^ - 2)gadgbc + (4(i + 10)fifacfi'f>d] + 2{d - 3)gacgbd 
^ [{d^ - 2)gabgcd + {d^ - 2)gadgbc + 2{3- d- d^)gacgbd\ 

^ ) [gabgcd + gacgbd + {d^ - "i) gadgbc] 



(A.26) 
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3o4 3o5 3o6 




Figure 6: The diagrams contributing to {ti{AiAj)tr{AkAi)) at order 
X'^/N'^, whose corresponding counterterms are relevant to the calculation of 



Mabcd = ( c^(rf+2) ) ^^^^ ~ ^)9acgbd + gab9cd + 9ad9bc\ 

^abcd^ = ^ 2) ) ^ ^)9ad9bc - C/ac^M - 9ab9cd] 

^ab^f = ^ ^ d{d + 2) ) ^ ~ 3ab9cd - 9ad9bc] (A.27) 

(3i) 5 / (i — 1 \ |. , 

^a6cd = \ d{d+2) J ^^"^^""^ ^ 9ac9bd + 9ad9bc\ 

^abcd = ( 4(i + 2) ) f^""^''^ ^ ^"'^^^^ ^ ~ 'i)9ab9cd] 



^abtd = ^ ( rf(rf + 2) ) - - fadi^bc 

(3Z) 5 / (i — 1 \ |. , 

^abcd = ^ \^ c^((i + 2) ) ''^"^^'"^ 9ac9bd + 9ad9bc\ 
^ibld ^ = ^ (^ d{d + 2) ) ^)9ac9bd - 9ab9cd - 9ad9bc 



^ilTd^^ - ^ (^ d{d + 2) ) ^^'^ ^)9ab9cd - 9ac9bd - 9ad9bc 



(A.28) 



^abcd = ^ (yd(JTij) ^^""^^^^ ^ ^"^^^'^ ^)9ab9cd] 
^abcd = ^ d{d + 2) ) ^^"^^^'^ + ^"^^^'^ -{d + l)9ad9bc\ 

^abcd = ( 373— ^ ) [9ab9cd + 9ac9bd - {d + l)9ad9b 



7r2 \d{d + 2)J (A.29) 

(3n4) _ f d- 1 

~ 47r2 I d 



Mabcd = — I ] 9ab9cd 



,^(3n5) 1 / - 1 1 



^'^abcd ~ ^'^abcd ~ ^^^abcd ~ ^^^abcd ~ ^^^abcd ~ ^^^abcd ~ ^■ 



We now sum over the appropriate permutations of indices to obtain the coefficient, 
R^jli, of the logarithmic divergence (|A.24| ) due to diagrams of each type. We hst these 
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coefficients, as well as the contribution to the tT{AiAj)tT{AkAi) counterterm i?^^ : 
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97r2 



(3n5) 
CT 



67r2 



In 



+ 



1 



367r^ 



(A.49) 



(A.50) 



Summing these yields the following result for the order X^/N"^ double-trace counter- 
term which is needed in equation (|3.48| ) : 



C 



CT 



In 



+ 



607r2 



[tr(A,A,)tr(AjAj) + 2tr(A,Aj)tr(A,Aj)] . 

(A.51) 



Appendix B. Useful Spherical Harmonic Identities 

In this appendix, we collect various properties of the 5''^ spherical harmonics required 
to study field theory on S"^ . Many of the basic results were derived in 0. 

B.l. Basic properties of spherical harmonics 

Scalar functions on the sphere may be expanded in a complete set of spherical har- 
monics SJ^'^ transforming in the (j/2,j/2) representation of SU{2) x SU{2) = SO{4), 
where j is any nonnegative integer, and — j/2 < m,m' < j /2. It is convenient to denote 
the full set of indices (j, m, m') by a. These obey an orthonormality condition (we take 
the radius of the 5""^ to be one) 

[ S'^S'^ = 5"^^, (B.l) 

JS3 

where 5"" denotes the complex conjugate of S"", 

^ gm m' ^ * ^ -j^^m+m' g—m —m' 2^ 

The spherical harmonics are eigenfunctions of the Laplace operator on the sphere, 

V25" = -ia(ia + 2)5«, (B.3) 
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and under a parity operation transform with eigenvalue (—1)^°'. 

A general vector field on the sphere may be expanded as a combination of gradients 
of the scalar spherical harmonics plus a set of vector spherical harmonics VJ!^. . These 
transform in the (^^^5 '^^) representation of 5'0(4), where j is a positive integer. Again, 
it is convenient to denote the full set of indices (j, m, m', e) by a single index a. These 
obey orthonormality relations 

(B.4) 

y« . V^^ = . 

S3 

Again indicates the complex conjugate of V", given by 

The vector spherical harmonics are eigenfunctions of parity with eigenvalue (—1)-'"'"^, and 
satisfy 

W" = -(i« + i)V", 

Vxy° = -e,(j, + l)y", (B.6) 
V ■ = 0. 

Explicit expressions for the scalar and vector spherical harmonics may be found in . 

B.2. Spherical harmonic integrals 

In expanding the action in modes, we require the integrals over the sphere of products 
of various numbers of spherical harmonics. For two spherical harmonics, the integrals 
are given by the orthonormality relations. For products of three spherical harmonics, we 
require the set of integrals given in ( |2.18| ). These were calculated in |^, and the results 
may be expressed in terms of the functions ( p.l2|) and (|3.13|) asH 



(jaP-f = ( 2 2 2(2 



ma mp rrij J V"^a ™ — ' 



J/3-e/3 




2 




m'^ 




J Ck. 


j 13-^13 


2 


2 







^"^^= 2 2 2 2 2 2 U3..e,(ja,J/3,i-y), 

^ / ia+e« j/3+g/3 j~f+<^-f \ f ja-Eg jp-<^P j',-<^-f \ 

^"^^= 2 2 2 2 2 2 U4e„e,e,0a,J/3,i7)- 

\ ttIq mp ruj J \ m J ' 

(B.7) 

The expression for C below differs by a factor of two from the expression in , but we beUeve 
that this expression is correct. 
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To evaluate integrals appearing in quartic terms, we do not require any additional infor- 
mation, since any product of two spherical harmonics may be expressed as a sum of single 
spherical harmonics using the completeness property and the integrals above. For example, 
we find 

V"" -V^ = D'^^^S^. (B.8) 

B.3. Identities involving 3j symbols 

The expressions for two and three loop vacuum diagrams involve products of the 
integrals in the previous subsection, with indices contracted and summed over in various 
ways. Since the m and m' indices appear only in 3j-symbols, we can always evaluate the 
sums over these using standard 3j-symbol identities. 

For basic manipulations, we require the identities 

jl h j3 \ ^ f 33 jl 32 \ /g 

mi m2 TO3 ) \TO3 mi 1712 J 

f jl 32 33 \ ^ ._-^x(ji+j2+j3) { jl j2 j3 \ /g 

\mi 1712 ms J \-mi -m2 -ma J 

f jl j2 j3 \ ^ /_-^^(Jl+J2+J3) f j2 jl js \ /g 

\mi m2 mz ) \m2 mi ms J 

To evaluate two-loop sums, we require 

( 31 j2 j3 \ f jl j2 jA \ ^ SjsjJmam^ /g -^2) 

^ ^ V mi 777,2 "73/1 "^1 "^2 "^4 / + 1 

mi m2 ^ ' ^ / ^ <j 

Finally, in three-loop computations we require 

^ ^^ ^_-^^Tni+Tn2+m3+m4,+Tn5+m6 

's 

jl jl is \ f is J4 J5 \ f is j4 je \ f J6 J2 J2 

mi —mi ms J \ —ms m4 m^ J \ — ms —7774 me / V ~™6 ^2 — m2 



= (-l)^-^+^V(2jl + 1)(2J2 + l)5,3,05,e,05. 
^ "^ ^_2^wi+"'i'2+m3+m4+m5+m6 



(B.13) 



Jl Jl -^2 W J2 J3 -^4 W J3 35 -^6 W J4 J6 J5 

mi — m2 m2 / V — m2 m3 m4 / I — ms m^ m^ ) \ — m4 —me — ms 



(B.14) 



57 



+m2 +m3 +m4 +7715 +7716 



31 J2 33 \ I J3 32 -^4 W J4 J5 -^6 W J6 J5 Jl 

mi 7712 W3 y y — ^3 —m2 1114^ J y — 7^4 7715 7775 J y — mg —7775 — 7r7i 



2ji + 1 



;B.15) 



mi +7772 +7773 +7774+7775 +7716 



Jl J2 -^3 W J3 J4 W J5 J6 ) ( -^2 J6 J4 

7771 W2 — TO3 / \?T73 7774 " ?T^5 / \ W5 7776 — TOl / \ — ?T^2 —7^4 

= C_l')Ji+i2+j3+j4+j5+j6|^l ^2 J3 I 

i4 is J6 

(B.16) 

In equations ( p.l4| ) and ( p.l5| ), the delta function with three arguments is either 1 or 



0, depending on whether or not the triangle relation is satisfied. In equation (|B.16|) , 

{ ^} } is the 6j symbol. 

3a 3b 36 



B.4- Identities for sums of spherical harmonics 

Using the 3j-identities, it is straightforward to derive expressions for sums over m,m' , 
and e in various products of the spherical harmonic integrals. For the two loop diagrams, 
we require: 

E D'^^'^D^^^ = 7^3c.{3c. + 2)jpU0 + 2)5-,,o, 
Yl D^^^D^P^ = ^Uj^ + 2)jp{jp + 2), 

m' s,e' s,j^ 

E D^^'^D^^^ = 2Rl^{j^,j,,jp) + 2Rl_{j^,j^,jp), (B.17) 

77i's,e's 

E D^f-^D^'P^ = ^Ja(ja + + 2), 

m' S,e' S,jry 

m' s 

For the three-loop diagrams it is useful first to note the basic relations (related to 
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those above): 



m's,e 


= —6 


XfijaUa + 2), 






m's,e's 


= 26yf 








m's,e's 






; Jr, ^7) + -R3- 


- {ja -1 jr ; J7 ) ) ; 


m's,€'s 










m's 




ec O'^: J7: Ja); 




m's 


= 0, 










= 0. 








m's 











(B.18) 



In each of these, we are summing over m, m' , and e for each of the contracted indices only. 

Using these basic relations and the results of the previous subsection, we find the 
following results for the non-vanishing spherical harmonic sums that appear in the three 
loop calculations: 

J2 i^-^A^7^A^75>^/3/3r ^ 5^.^^^.^ ^^^^j^ ^ 2)j^(j^ + 2)j,{j, + 2), 



m's,e's 



■ {Rs+{jj. jrjp) + R3- UjJrJfi)), 

J2 DC.,Xj^^,rj^5p-Xj^-5-pf ^ 5j,,,^—^{Rl+{j^,jx,j,)+RLUa,jx,j,)) 
n's,e's ^-^^ ^ 

• (Rl+UiaJxJs) + Rl-Ui3,jx,js)), 

\^ jjaPXjjjSXjjPSr jja^yr _ f-^\2_^j s \ 2 2 2IJ2 2 2I 



Rsec^e^ (jajx, 3p)Rze^es {j^i : jx, 35)Rz€pes UpJt, js)R3ece-y (jajri 

(B.19) 
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e's 



4 



J2 C^^PC-P-^^D-P^D^^- = -5,^,j^j-—r^Rl{3o.3,,3p) 

n's,e's ^ ^> 

■ iR3+{ja,ja,j(3) + Rs-iJa, ja, jp)), 

7i's,e's '^^ ^ 

■^{R4epesep{j(3,js,jp)f, 
e's 



e~f e's 



(B.20) 

ja+Eg Jg+e/3 J-r+£T ^ ( ja-<^a 3l3—<^l3 3r—<^-, 



„ , r Ja+eg Jg+e/3 J-t+Et ^ r Jg-Ea Jp-<^P 3t-<^t ~\ 

EnajPnPSpjpa/STTpS'yT _ _(_-,\2_^j s I 2 2 2IJ2 2 2I 

1 J,5+£6 i^ + ^T 2£. I 1 35-^6 3i-^i 3p_ { 

J ^ 2 2 2-* 



■ Rze^e^ {.joii jpi jy)R3e0es U/Si jpi jd^RAeaepeT {,jai j Pi jr^RAese^eT {jSi j^i jr^ 

^_ - _ / r -^P+^P JcT+eg jT+<^T "I ( 3p-^p icr-ec Jt-^t ^ 

T^paT t^ocPt T^J3'ipT^-ia(T _ s \ 2 2 2lJ2 2 2l 

±J ±J ±J ±J y ±J N j^^g^ j/3+e/3 ^7+^7 I I 3a-ea 3/3-^0 j-j-^l ( 

m's '-2 2 2>''-2 2 2>' 

• R^epe^e^ Up 1 jcn jr^RAeaepeT {joij j/St j t) R^e fs e^ e p {jpi J75 j p^RAe^e^e^ U^t 3a-i 3a)-i 

^ _ _- , r 2£. k. 313+^13 ^ f 2p It 3fi-<^i3 ^ 

■Qa-yp-Q^'yaQpPT^Taa _ _/-^\2_^j s \ 2 2 2IJ2 2 

/ ^ ^ 1 3~1+'^~1 ja+£a I 1 is: 3^-^^ Jet -Eg ( 

m's ^ 2 2 2>''-2 2 2'' 

• -R2(jp, J>, jr)R2{jTjc.Ja)R3e^e^ {jcjp, j'r)R3epe^ Ja, J7), 

^ _ __ ^ f jg+eg jp+<^() 3-1+^1 •\ ( jg-ea 3f3-<^f3 3-,-^! ^ 

E^aP-y/jpja jj^ra jjarp _ (_l)2_^3 s \ 2 2 2IJ2 2 

I is: ip. 3T+e-T I I is: ip ir-er | 

m's ^-2 2 2>''-2 2 2>' 

■ R2{jpi j'yt ja)R4eaepe^ {jai jf3i jl^R^epe-r {jpi jaj jr^R^eaeT {jai jpi 3t)- 

(B.21) 

Finally, in certain cases, we may simplify expressions further by using the relations 

Rl+{a, c,b) = J2 Rl- c, h) = -^a{a + 2)6(6 + 2), 

(B.22) 

^(i?2(a, b, c)f = ^b{b + 2)a{a + 2)(a + 1)2. 
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